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ACOUSTIC  INTENSITY  VECTORS  FROM  AN  INFINITE  PLATE 
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-^Formulae  are  given  for  the  acoustic  pressure  and  particle  velocities  from 
which  the  acoustic  intensity  vectors  are  calculated.  Plots  of  the  vectors 
illustrate  their  utility  as  a  visual  aid  to  the  understanding  of  the  physics 
of  fluid-structure  interaction. 
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LIST  OF  SYMBOLS 


(x,z) 

Cartesian  coordinates 

Wfx) 

plate  displacement 

W(x,z),  U(x,z) 

fluid  particle  displacements 

W(x,z) ,  U(x,z) 

fluid  particle  velocities 

p(x,z),  p*(x,z} 

acoustic  pressure  and  its  complex  conjugate 

p  (o,z) 

Fourier  transforms  of  plate  displacement  and 

D 

fl3xural  rigidity  of  plate,  «Eh3/12(1-v2) 

E,  n 

Young's  modulus  and  loss  factor 

V 

Poisson’s  ratio 

V  h 

density  and  thickness  of  plate 

P,  C 

density  and  sound  velocity  of  fluid 

f 

frequency  of  line  force,  Hz 

u 

radian  frequency,  «2*f 

k 

acoustic  wavenumber,  *2nf/c 

Fo 

magnitude  of  line  force 

1 

Dirac  delta  function 

Ix(x,z,t),  Iz(x,z,t) 

instantaneous  acoustic  intensities 

yx,z),  I2(x,z) 

time  averaged  acoustic  intensities 

M,  S 

mass  and  spring  constant 

{W} 

n^1  displacement  matrix 

(F),  {P} 

nx1  "x1 

n  1  force  matrices 
* 

ID],  (D]~l 

nxn  nxn 

n^n  receptance  matrix  and  its  inverse 

[A] 

nxn 

n  n  diagonal  matrix  with  elements  a.. 
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INTRODUCTION 


The  time  averaged  acooatio  intensity  vector  diagrams  not  only  provide 
concise  visual  information  about  the  magnitude  and  direction  of  the  energy  flow 
in  an  acoustic  fluid*  but  also  serve  as  an  aid  to  the  understanding  of  the 
complex  physios  involved  in  fluid-structure  interactions.  Acoustic  intensity 
vectors  oan  be  determined  either  experimentally  or  theoretically.  Fahy  (1]  has 
used  a  prototype  intensity  meter  to  measure  the  intensity  distribution  around  a 
small  machine,  thus  demonstrating  the  practicality  of  measurement.  Other 
authors  [2.3]  have  suggested  procedures  for  measuring  the  acoustic  intensity 
close  to  a  vibrating  surface.  Williams  et  al  [4]  have  constructed  a  microphone 
array  of  256  elements,  whioh  is  used  in  the  nearfield  of  a  vibretori  subsequent 
computer  processing  enables  reconstruction  of  the  acoustic  field  and  particle 
velooities.  and  hence  acoustic  Intensity  vectors.  Kristiansen  [5]  has  produced 
theoretical  vectors  due  to  assumed  mode  shapes  of  a  membrane  vibrating  in  an 
infinite  baffle. 

Here,  acoustic  intensity  vectors  from  an  elastic  plate  with  line  attach¬ 
ments  are  derived  theoretically.  Thr;  line  attachments  may  be  time-harmonic 
forces,  masses,  masses  and  springs,  ground°d  springs  or  displacement  constraints. 
In  Section  2  formulae  are  given  for  the  acoustic  pressure  and  particle  veloci¬ 
ties  due  to  line  force  excitation  of  the  plate.  The  acoustic  intensity  vector 
is  defined  in  Section  3.  The  procedure  necessary  to  calculate  the  dynamic 
stiffness  of  the  plate  with  line  attachments  is  given  in  Section  4,  where  the 
equivalent  forces  due  to  the  line  attachments  are  also  calculated.  Finally,  in 
Section  5,  acoustic  intensity  vector  plots  for  forces  and  constraints  are 
discussed. 


2.  ACOUSTIC  PRESSURE  AND  DISPLACEMENT  DUE  TO  LINE  FORCE 

A  uniform  infinite  thin  plate  occupies  the  plane  z=0,  and  it  is  subjected 
to  a  prescribed  time-harmonic  line  force,  FQexp(-iwt),  located  at  x=xo.  One 

side  of  the  plate  is  in  contact  with  an  acoustic  fluid  filling  the  half-space 
z>0,  and  the  other  side  contacts  a  vacuum.  The  geometry  of  interest  is  shown 
in  Figure  1.  Time  variation,  exp(-iwt),  will  be  omitted  from  all  equations. 

The  plate  displacement  and  fluid  sound  pressure  are  expressed  as  Fourier 
transforms  t 


Wfx)  •  (1/2ir) 


p(x,z)  B  (1/2tt) 


2ir)  I  W(o 

i  —CO 

2ir)|"  'pCo 


)exp(ictxjda 


,z)exp(iax)da 


Substitute  these  transforms  into  the  plate  and  accustic  equations  of  motion 

[  D34/3x4-w2p  h]W(x)  =  F  $(x-x  )-p(x,0)  C 

5  u  U 


(32/3x2  +  32/3z2)p(x,i:)  =  -k2p(x,z) 


to  give,  after  carrying  out  some  straightforward  algsbra. 


(S) 


W(x)  ■  (F„/2w)  ■»Pll°(»-‘<o))da 

h-ipw^/Y 

8 

p(x,z)  -  C-WF/2.)  f  «xp»rri-xpll=.(x-x0Ud. 

**•  Y(Da4-d)2p  K-ipw2/y) 

8 

where  y  *  Ak2-a2)  with  Im(Y)>0  in  order  to  satisfy  the  radiation  condition. 
It  is  to  be  noted  that  the  normal  fluid  particle  displacement  (cf  Eq.  7)  at  the 
plate  surface  is  equal  to  the  plate  displacement.  This  is  the  boundary 
condition  of  normal  displacement  continuity. 


The  acoustic  particle  displacements  in  the  fluid  are  given  by  the  equations 
pw2  lU(x,z),W(x,z)J  ■  l3p(x,z)/3x,3p(x,z)/3z]  (7) 

which  beoomB,  on  making  use  of  Eq.  6. 


U(x,z) 


W(x,z) 


(F  /2ft)  a«B*P(iVz)exptitt(x-xo))da 
0  i_  Y(Da4-w^P  h-ipoi2/Y) 

(F  /2t r)  exp(iYz)exp[ia(x»x0)]da 

0  [_  Da4-u2p  h-ipw2/Y 

*'*'  S 


(0) 


The  acoustic  pressure  and  fluid  r  rticle  displacements  due  to  a  number  of 
line  forces  cun  be  found  by  linear  superposition. 

3.  ACOUSTIC  INTENSITY  VECTORS 


Let  p(x,z),  U(x,z)  and  W(x,z)  be  respectively  the  complex  pressure, 
fluid  particle  x-displacsment  and  fluid  particle  z-displacement  due  to 
externally  applied  line  forces.  The  fluid  particle  velocities  are  derived  by 
time-differentiating  the  displacements.  They  are 

U(x,z)  ■  -iu>U(x,z) 

W(x,z)  ■  -iwW(x,z)  (9) 

The  instantaneous  acoustic  intensities  are  defined  as 

e 

Ix(x,z,t)  ■  Re  tp(x,z)exp(-iwt)]  .Re  (U(x,z)exp(-iu)t)] 

Iz(x,z,t)  ■  Re  tp(x,z)exp(-iu)t)i  .Re  [W(x,z)exp(-iut;] 


MO) 


The  time  averaged  acoustic  intensities  are  given  by 


r  T 

T  (x,z)  ■  limi 
* 

Ix(x,z,t)dt 

T-*-  J 

0 

T  (x,z)  ■  lim^ 

^  T 

rT 

Iz(x,z,t)dt 

T-h» 

0 

which  reduce  to  the  formulae 

e 

Ix(x,z)  -  iRetp*(x,z)U(x,z)l 

e 

I2(x,z)  -  iRe tp*(x,z)W(x,z)l 


(11) 


(12) 


The  amplitude  and  direction  of  the  resultant  acoustic  intensity  vector  are 
respectively 


I(x,z)  -  /[Ix2(x,z)  +  Iz2(x,z)] 

0  ■  tan“*(I  (x,z)/I  (x,z))  (13) 

4  A 

4.  EQUIVALENT  FORCES  DUE  TO  LINE  ATTACHMENTS 


(a)  General 


The  plate  dynamic  stiffness  equation  must  first  be  formed.  From  this, 
the  system  dynamic  stiffness  equation  is  then  derived  using  standard  finitB 
element  methods.  This  system  equation  is  solved  to  find  the  plate  dis¬ 
placements,  which  than  yield  the  equivalent  forces. 

(b)  Plate  Dynamic  Stiffness  Equation 

Let  n  separate  line  forces,  (F),  act  on  the  plate  at  the  points 
Xj,  i"1,n.  Then  the  n  plate  displacements,  (W),  are  related  to  these 
forces  by  the  frequency  dependent  receptance  matrix,  [D] ,  where 


(W)  -  [D]  .  {F} 


(14) 


"x1 


n  n  n  1 
x  x 


The  ij  "  element  of  [D]  is  simply  the  response  at  x^  due  to  unit  force 
at  x..  It  is  given  by  Eq.  5.  The  inverse  of  Eq.  14  is  the  plate  dynamic 

w 

stiffness  equation. 


ID]-1  .  {W}  =  {F>  (15) 

n  n  n  1  n  1 
x  xx 


S 
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(o)  System  Dynamic  Stiffness  Equation 

The  system  dynamio  stiffness  squation  is  of  the  form 

lD_l  +  A]  .  {W}  -  (P>  (16) 

V  "x1  "x1 

where  [A]  is  a  diagonal  matrix  whose  elements  aro  derived  from  the  mass/ 
spring  attachments,  and  (P>  is  a  oolumn  matrix  of  prescribed  line  foroes. 
The  mass/spring  attachments  are  incorporated  as  follows i 

For  a  mass  located  at  x^,  set  a^  to  -w2M. 

For  a  grounded  spring  at  x^,  set  to  S. 

For  a  mass  and  spring  at  x^,  set  a^  to  -w2M/(1-w2M/S) . 

A  displacement  constraint  at  x,  is  represented  by  setting  the  ith  row 

A.L.  1 

and  itn  column  of  [D] “1  to  zero,  and  a^  to  unity. 

(d)  Equivalent  Forces 

The  system  dynamic  stiffness  equation  is  inverted  to  give  the  plate 
displacements,  {W>,  which  are  inserted  into  the  plate  dynamic  stiffness 
equation,  Eq.  15,  to  yield  the  equivalent  forces,  {F}.  The  acoustic 
pressure  and  fluid  particle  velocities  can  then  be  found,  as  described  in 
Section  2. 

NUMERICAL  RESULTS 


(a)  General 


Two  Fortran  programs  have  been  written  which  together  computs  and  plot 
the  intensity  vectors.  The  first  program  calculates  and  stores  on  a  disk 
file  the  pressure  and  particle  velocities,  at  points  on  a  selected  grid  in 
the  positive  (x,z)  quadrant,  duo  to  a  single  line  force  excitation.  The 
second  program  uses  this  information  to  compute  and  plot  the  acoustic 
intensity  vectors  due  to  the  plate  with  line  attachments.  Examples  of  the 
graphical  output  are  shown  in  Figures  2-7  for  which  the  following  values  in 
SI  units  were  used: 

plate:  E  »  19.5E10,  v  ■  0.29,  p  *  7700.0,  h  -  0.02 

s 

fluid:  p  -  1000.0,  c  -  1500.0 

Plate  damping  is  represented  by  a  complex  Young's  modulus,  E(1-in), 
where  n,  the  loss  factor,  is  chosen  as  0.02,  The  plotted  vector  lengths 
are  proportional  to  I(x,z)  in  Figure  5,  and  to  /fTCx.z)]  in  the  other 
intensity  plots. 

(b)  Single  Line-Force  Excitation,  250  Hz 

Figures  2  and  3  show  the  intensity  plots  which  illustrate  the  essen¬ 
tial  acoustic  features  due  to  low  frequency  excitation.  First,  the  sound 
radiation  is  very  inefficient.  Secondly,  the  intense  near  field  shows 
areas  where  energy  is  bsing  received  by  the  plate.  Thirdly,  the  fluid 
surface  wave  associated  with  the  subsonic  plate  flexural  wave  is  evident. 
Finally,  a  closed  loop  of  energy  circulation  can  bs  seen. 


(o)  Single  Lina-Forca  Exoitation,  2.5  kHz 


Figure  4  shows  the  intensity  plot  at  a  mid-frequency  which  illustrates 
the  increasing  radiation  effioienoy  of  the  plate. 

(d)  Single  Line-Force  Excitation.  25  kHz 

Figure  5  shows  the  intensity  plot  at  a  high  frequency.  At  this 
fraquanoy,  whioh  is  approximately  twice  the  so-called  critical  frequency, 
the  'leaky*  wave  beam  [6]  can  be  seen  emanating  from  the  plate  surface. 
The  attenuation  along  the  beam  is  small. 

(e)  Effsot  of  Constraints.  250  H2 

Figures  6  and  7  show  the  vectors  due  to  a  force  located  between 
symmetrically  placed  displacement  constraints.  In  Figure  6  the  system  is 
at  its  fundamental  resonance i  the  top  row  of  intensity  vectors  shows 
radiation  coming  from  the  region  of  the  constraint.  In  Figure  7  the 
constraints  cause  intense  near-fields  with  apparently  little  radiation i 
it  is  interesting  to  nets  that  energy  enters  the  plate  at  the  drive  point. 

6.  ACKNOWLEDGEMENTS 


Thanks  are  due  to  J  H  Games  for  his  helpful  supervision  of  the  project, 
and  to  E  3  ClemBnt  for  his  advice  on  the  computer  plotting  facilities. 


W.J,  Spicer  (Student,  Mathematics  Department, 
College  of  St  Hild  and  St  Bede, 
Durham  University) 


WJS/jms 


REFERENCES 


FAHY,  F.J.,  Measurements  with  an  intensity  meter  of  the  acoustic  power  of 
a  small  machine  in  a  room,  0,  Sound  Vib.,  57(3),  April  1975,  pages  311-322. 

HANSEN,  C.H.,  BIES,  D.A.,  N'W'  field  determination  of  the  corrplex  radiation 
efficiency  and  acoustio  inte  ,sity  distribution  for  a  resonantly  vibrating 
surface,  3.  Sound  Vib.,  62(1),  January  1979,  pages  93-110. 

PETTERSEN,  O.K.O.,  A  procedure  for  determining  the  sound  intensity  distri¬ 
bution  close  to  a  vibrating  surface,  3.  Sound  Vib.,  66(4),  October  1979, 
pages  626-629. 

WILLIAMS,  E.G.,  MAYNARD,  3.D.,  SKUDRZYK,  E.,  Sound  source  reconstructions 
using  a  microphone  array,  J.Acoust.Soc.Am. ,  66(1),  July  1960,  pages 
340-344. 

KRISTIANSEN,  U.R.,  A  numerical  study  of  the  acoustic  intensity  distribution 
close  to  a  vibrating  membrane,  3. Sound  Vib.,  76(2),  June  1961,  pages 
305-309. 

CRIGHTON,  D.G.,  The  free  and  forced  waves  on  a  fluid-loaded  elastic  plate, 
3.  Sound  Vib,,  63(2),  March  1979,  pages  225-235. 


ORIVE 

FORCE 


MASS  ANO 
SPRING 


GROUNOEO 

SPRING 


DISPLACEMENT 

CONSTRAINT 


FIG  .  1 


PLATE  COORDINATE  SYSTEM  AND  ATTACHMENTS 


AMTC ( N) TMEIOEE 


-  13  ' 


FIB*  2  2S0Hx  DRIVE  FORCE  AT  X-0.0 
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FIG.  3  250Hz  DRIVE  FORCE  AT  X-0.0 
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FIG.  6  DISPLACEMENT  CONSTRAINTS  AT  E0,2»  FROM  DRIVE  FORCE.  250Hz 
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